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The effects of electron doping and phonon vibrations on the magnetic properties of monolayer 
and bilayer FeSe epitaxial films on SrTi 03 have been studied, respectively, using first-principles 
calculations with van der Waals correction. For monolayer FeSe epitaxial film, the combined effect 
of electron doping and phonon vibrations readily leads to magnetic frustration between the collinear 
antiferromagnetic state and the checkerboard antiferromagnetic Neel state. For bilayer FeSe epi¬ 
taxial film, such magnetic frustration is much more easily induced by electron doping in its bottom 
layer than its top layer. The underlying physics is that the doped electrons are accumulated at 
the interface between the FeSe layers and the substrate. These results are consistent with existing 
experimental studies. 

PACS numbers: 74.20.Pq, 74.70.Xa, 74.78.-w, 75.70.Ak 


I. INTRODUCTION 

Iron-based superconductors have drawn worldwide in¬ 
terests both experimentally and theoretically since their 
discoveries ji - — The superconducting transition tempera¬ 
ture T c can be modulated by charge doping via chem¬ 
ical substitution or by high pressure. In addition, epi¬ 
taxial film growth is another effective way to tune the 
lattice parameters, magnetic properties, and T c of Fe- 
based superconductors £ - — By using molecular beam epi¬ 
taxy (MBE) technique, atomically flat FeSe epitaxial 
films have been grown on bilayer graphene^ and SrTiOs 
substrates^ Surprisingly, the monolayer FeSe on SrTiOs 
shows a superconducting transition signature around 77 
K, exceeding the highest T c of bulk ferropnictides^ while 
the bilayer FeSe does not exhibit a superconducting gap 
in the scanning tunneling microscope (STM) spectral 
With larger tensile strains on FeSe epitaxial films, the 
T c can also be lifted to ^ 70 Due to its simplest 

atomic structure among Fe-based superconductors, these 
fresh experiments have demonstrated that FeSe epitaxial 
film is a good model system for exploring the unconven¬ 
tional superconducting mechanism in Fe-based supercon¬ 
ductors. 

Various experiments have been carried out to inves¬ 
tigate the physical properties of FeSe epitaxial films 
on SrTiOs ii - — Concerned with the electronic struc¬ 
ture of FeSe monolayer on SrTiOs, only electron-like 
Fermi surfaces exist at the corners of the Brillouin 
zone (BZ) in the angle resolved photoemission spec¬ 
tra (ARPES) measurement ii - — This is quite different 
from that of multilayer FeSe films and bulk FeSe. It 
has been suggested that charge transfer from the sub¬ 
strate to FeSe layers plays an important role in the 
superconductivity*^ - — For both monolayer and bilayer 
FeSe, more and more electrons are transferred from the 
substrate when annealing time increases.- 0 As for the 
magnetism, it was demonstrated that superconductivity 
occurs when the electron transfer from the substrate sup¬ 


presses the otherwise pronounced spin density waves in 
monolayer FeSeJ^ Recently, replica bands with a dis¬ 
tance of 100 meV in ARPES measurement were found 
for monolayer FeSe epitaxial filing It was then proposed 
that this is induced by the coupling between electrons 
and substrate phonons, which enhances the supercon¬ 
ducting pairing temperature profoundly*^ The above in¬ 
tensive experiments have provided valuable information 
about the FeSe epitaxial system. 

Many efforts have also been devoted from a theoretical 
standpoint*^ - — In our previous calculations^ the sim¬ 
ilar electronic structures of monolayer and bilayer FeSe 
on undoped SrTiOs suggested that the superconductiv¬ 
ity would occur at the first layer of FeSe epitaxial film or 
at the interface. Xiang et ai studied the influence of the 
screening effect of a ferroelectric phonon of the substrate 
on the Cooper pairs in the FeSe layer. 23 The electron dop¬ 
ing effect has also been investigated*^ - — The measured 
Fermi surface of FeSe monolayer can be reproduced in 
calculation when the antiferromagnetic (AFM) Neel or¬ 
der is set np j 24 i 25 Cao et ai have studied the effects of 
tensile strain and charge transfer on the spin density wave 
in FeSe/SrTiOs thin films*2& Li et al. have calculated the 
electron-phonon coupling constant of FeSe/SrTiOs using 
the first-principles approach, indicating that the electron- 
phonon mechanism alone cannot explain its high T C S^ 
The above theoretical works have proposed a variety of 
possible contributions to the superconductivity in FeSe 
epitaxial film from different points of view. Nevertheless, 
the physical mechanism of superconductivity in this sys¬ 
tem is still an open question. In particular, is the high 
transition temperature correlated with magnetism? Is 
there any other substantial phonon effect if the electron- 
phonon coupling alone cannot account for the high T c ? 
The entanglements of tensile lattice strain, electron dop¬ 
ing, phonon, and magnetism in the FeSe epitaxial film 
further complicate the situation. 

In this paper, we have studied the effects of electron 
doping and phonon vibrations on the magnetic properties 
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of FeSe monolayer and bilayer on SrTiOs by using first- 
principles calculations with van der Waals (vdW) cor¬ 
rection. The nonmagnetic, checkerboard AFM Neel, and 
collinear AFM states of FeSe epitaxial films have been in¬ 
vestigated. The relationships of AFM fluctuations with 
electron doping and quantum zero-point (ZP) atomic dis¬ 
placements of phonon modes have been examined. The 
differences between monolayer and bilayer FeSe epitaxial 
films have also been addressed. 

The rest of this paper is organized as follows. In Sec. 
II, the computational details are described. In Sec. Ill, 
the AFM variations of FeSe epitaxial films on SrTiOs 
have been studied as functions of electron doping and 
ZP atomic displacements. Discussions with related ex¬ 
perimental and theoretical works are given in Sec. IV, 
and a short summary is provided in Sec. V. 


II. COMPUTATIONAL DETAILS 

Fully spin-polarized first-principles calculations were 
carried out by using the projector augmented wave 
(PAW) method^ as implemented in the Vienna ab- 
initio simulation packaged The generalized gradient ap¬ 
proximation (GGA) of Perdew-Burke-Ernzerho® for the 
exchange-correlation potentials was adopted. To describe 
the vdW interaction in layered systems not included in 
the conventional density functional theory, the vdW- 
optB86b functional^ was chosen. To model FeSe ultra- 
thin films on Ti02-terminated SrTiOs (001) as grown in 
experiment,- we used a six-layer SrTiOs(001) slab with 
FeSe monolayer and bilayer adsorbed on the top side in a 
a/ 2 x \J~2 two-dimensional supercell plus a vacuum layer 
> 10 A. The kinetic energy cutoff of the plane-wave ba¬ 
sis was chosen to be 400 eV. An 8x8x1 k-point mesh for 
the Brillouin zone sampling and the Gaussian smearing 
technique with a width of 0.05 eV were used. The non¬ 
magnetic, checkerboard AFM Neel, and collinear AFM 
states were studied. In all these magnetic orders, the top 
two slab layers and all atoms in FeSe layer(s) were al¬ 
lowed to relax until the corresponding forces were smaller 
than 0.01 eV/A, while the bottom slab layers were fixed 
at their bulk positions. The electric field induced by 
asymmetric slab relaxation was compensated by a dipole 
correction^ 

We note that as FeSe has charge-neutral Se-Fe 2 -Se lay¬ 
ers, the vdW interaction plays an important role in the 
interlayer bonding. Our previous studies on bulk FeSe 
have shown that the lattice parameters (especially along 
the c direction) and phonon frequencies can be accurately 
calculated only when the vdW interaction is taken into 
account^ In the present study on FeSe epitaxial system, 
we had checked the equilibrium distance between the top 
Se atom of the first FeSe layer and the Ti02 termination 
layer of the substrate. The calculated distance with vdW 
correction (5.56 A) is in excellent agreement with the ex¬ 
perimental value (5.5 A) #2 In comparison, the calculated 
distance without vdW correction is 5.85 A. The inclusion 


of vdW correction in the calculations is very important 
to get realistic adsorption structural parameters, which 
is a prerequisite to study the related properties such as 
the phonon and its derivatives. 

After the equilibrium structures were obtained, the 
frequencies and displacement patterns of phonon modes 
were calculated using the dynamical matrix method. 37 
Firstly, we calculated the phonon modes with equilibrium 
structure in the collinear AFM order (the energetically 
favorable state). Secondly, even though the ground state 
is in collinear AFM order, magnetic domains will always 
form in a real material with magnetic stripe directions 
randomly along the a or b direction. It turns out that 
the ARPES measured Fermi surface is an average effect 
of the electronic states of those magnetic domains. Cor¬ 
respondingly, in calculations, the Fermi surface in the 
checkerboard AFM Neel state can reproduce the ARPES 
results ] 24 ! 25 Thus we also studied the lattice dynamics in 
the checkerboard AFM Neel order. Our previous study 
on bulk FeSe has shown that the phonon frequencies cal¬ 
culated using the checkerboard AFM Neel order fit the 
experimental results well^ We have checked that the dif¬ 
ferent antiferromagnetic orders do not change the phonon 
frequencies very much (< 2 meV). In total, there were 
54 and 78 independent phonon modes calculated for the 
monolayer and bilayer FeSe epitaxial films, respectively. 

Physically the influence of a phonon on the magnetic 
properties of FeSe epitaxial film happens through the 
electron-phonon coupling, which has twofold impacts on 
electronic structure, namely thermally scattering of the 
electron and quantum zero-point vibration of the phonon. 
Here we study the impact resulting from the zero-point 
vibrations of phonon modes. In the calculations, the 
atomic displacements due to the zero-point vibrations of 
phonon modes were obtained according to the method of 
Ref. 0. Specifically, the atoms were displaced to a vibra¬ 
tional state with a potential energy of hu s /4 instead of 
Tiujs/2 for each specified phonon mode s according to the 
Heisenberg uncertainty principle, while its normal-mode 
coordinates could reach two maxima along two opposite 
directions^ When the atoms were displaced, the ener¬ 
gies in different AFM orders were calculated on the same 
distorted structure. The variation of the local magnetic 
moment on Fe due to zero-point atomic displacement is 
| A M\ = |M_|_ — M- 1, with M + being the local magnetic 
moment on Fe for one displacement and M_ the other 

Here the electron doping effect can be simulated with 
either oxygen vacancies in the substrat o 26 ! 28 or excess 
electrons in the system. Considering that the oxygen va¬ 
cancies may induce some artificial vibrational modes, we 
thus simulated the electron doping effect by changing the 
total number of electrons in the system, which reserves 
the spatial symmetry as the undoped case. As to the 
amount of doped electrons, the ARPES experiment sug¬ 
gests an electron counting of about 0.1 electrons/Fe by 
analyzing the area ratio of the Fermi surface around M to 
the whole Brillouin zone*22 Accordingly, we studied the 
systems with electron doping from 0.1 to 0.4 electrons for 
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TABLE I. The relaxed layer distances d (in A) and local magnetic moment M (in fi b) on Fe of monolayer and bilayer FeSe 
epitaxial films in the collinear AFM order and the checkerboard AFM Neel order under different levels of electron doping. The 
Mi and M 2 denote the magnetic moment of Fe in the bottom and top FeSe layers, respectively. The ’Neel-Collinear’ means 
the bottom layer of FeSe bilayer is in the checkerboard AFM Neel order while the top layer in the collinear AFM order. 
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FIG. 1. (Color online) Atomic structures of (a) monolayer 
and (b) bilayer FeSe on TiO 2 -terminated SrTiOs (001) sur¬ 
face. Magnetic patterns of Fe atom spins denoted by red ar¬ 
rows in the (c) checkerboard AFM Neel state and (d) collinear 
AFM state. 


the whole supercell. The charge difference density was 
calculated from A Q = Q(FeSe/SrTiOs) - Q(SrTiOs) - 
Q(FeSe), where Q(SrTiOs) is the charge density of the 
undoped or electron-doped substrate. 


III. RESULTS AND ANALYSIS 

The structures and spin patterns of monolayer and bi¬ 
layer FeSe films on SrTiOs are shown in Figure [lj As 
in the MBE experiment^ FeSe films are adsorbed on 
the Ti02-terminated SrTiOs (001) surface. From the to¬ 
tal energy calculations, we find that the FeSe epitaxial 
film energetically favors such adsorption sites at which 
the bottom Se atoms are on top of the terminated Ti 
atoms of the substrate [Fig. QJa) and QJb)]. This is 
also consistent with our previous calculation^ For mag¬ 


FIG. 2. (Color online) Energy differences between differ¬ 
ent magnetic states for monolayer FeSe on Ti02-terminated 
SrTiO3(001) surface as functions of electron doping. 


netic orders, the nonmagnetic, ferromagnetic, checker¬ 
board AFM Neel, and collinear AFM states have been 
considered. The spin patterns of the latter two cases are 
shown in Fig. QJc) and[]Jd), respectively. The relaxed 
layer distances d and local magnetic moments M on Fe 
in the collinear AFM order (the energetically favorable 
state) and the checkerboard AFM Neel order under dif¬ 
ferent levels of electron doping are listed in Table I. For 
comparison, the structural parameters of the FeSe epitax¬ 
ial films in the nonmagnetic state are provided in Table 
II in the appendix. 


A. Monolayer FeSe film on SrTiOs 

Figure [2] shows the calculated energy differences be¬ 
tween checkerboard AFM Neel and nonmagnetic states 
as well as between collinear AFM and checkerboard AFM 
Neel states as a function of electron doping. As we 
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FIG. 3. (Color online) (a) Three dimensional and (b) one 
dimensional charge difference density for monolayer FeSe ad¬ 
sorbed on Ti02-terminated SrTiO3(001) surface with 0.2- 
electron doping. Dark red and blue isosurfaces in panel (a) 
are respectively electron accumulation and depletion areas in 
isovalue of 0.01 e/A 3 . 


can see, the checkerboard AFM Neel state remains en¬ 
ergetically lower than the nonmagnetic state in all the 
doped cases, and the corresponding energy difference 
only changes slightly, within 2 meV/Fe. In contrast, the 
energy difference between the collinear AFM state and 
the checkerboard AFM Neel state decreases considerably 
when electrons are doped. To be specific, the energy of 
the collinear AFM state is 23.5 meV/Fe lower than that 
of the checkerboard AFM Neel state in the undoped case, 


while it decreases to 17.4 meV/Fe in the case of 0.2 elec¬ 
trons doping. However, after more than 0.2 electrons are 
doped, the energy difference no longer decreases. Thus 
the electron doping indeed has a notable effect on the 
magnetic properties of monolayer FeSe epitaxial film, but 
it is still not enough in and of itself to change the mag¬ 
netic states energetically. 

To clarify how the charge transfer takes place between 
the doped substrate and the FeSe adlayer, we inspect 
the calculated charge difference density in the case of 
0.2-electron doping, as shown in Fig. [3] Clearly, the 
doped electrons converge to the interface between the 
FeSe adlayer and the SrTiOs substrate. Meanwhile, the 
adsorption of the FeSe layer on the substrate induces 
a certain amount of charge redistributed around the Fe 
atoms, which is responsible for the energy changes among 
the different AFM states. Interestingly, a recent study 
has found strong coupling of the Iron-quadrupole and 
anion-dipole polarizations in Ba(Fei_ a ,Co ;E )As 2 . 39 Phys¬ 
ically the electron accumulation at the interface provides 
the bonding interaction between the adlayer and the sub¬ 
strate, which is expected to be very susceptible to the 
inter facial atomic vibrations. 

In our previous study on crystal /3 -FeSe under hydro¬ 
static pressure^ it was found that in comparison with 
the other phonon modes, the zero-point atomic displace¬ 
ment due to the coherent vertical vibration of Se induces 
the largest variation of local magnetic moment on Fe4& 
This urges us to examine whether or not the same phonon 
effect exists in FeSe epitaxial films. Figures 0] and [5] show 
the calculated variations of the local magnetic moment in 



0 5 10 15 20 25 30 35 40 45 50 55 

Label of vibrational modes 


(a) 


undoped 


0.20 
0 . 15 - 
0 . 10 - 
0 . 05 - 

3 

Q.oo ltttt , , as , gtt 


□ |AM|(Col) 
A |AM| (Neel) 


□ 

A 


a A ° „ 

a ■ A A Afo, A 



Label of vibrational modes 


FIG. 4. (Color online) Variations of local magnetic mo¬ 
ment on Fe induced respectively by the zero-point atomic 
displacements of all the vertical phonon modes for (a) un¬ 
doped and (b) 0.2-electron doped monolayer FeSe on SrTi 03 . 
The phonon modes calculated at the equilibrium structure of 
the collinear AFM state. 


FIG. 5. (Color online) Variations of local magnetic mo¬ 
ment on Fe induced respectively by the zero-point atomic 
displacements of all the vertical phonon modes for (a) un¬ 
doped and (b) 0.2-electron doped monolayer FeSe on SrTiOs. 
The phonon modes calculated at the equilibrium structure of 
the checkerboard AFM Neel state. 
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FIG. 6. (Color online) Energy difference between collinear 
AFM and checkerboard Neel AFM states varying with the 
positive (p) or negative (n) zero-point atomic displacements 
of vertical phonon modes for monolayer FeSe on SrTi 03 in (a) 
undoped and (b) 0.2-electron doped cases. Grey areas high¬ 
light the regions with energy difference within ±10.0 meV/Fe. 
The phonon modes calculated at the equilibrium structure of 
the collinear AFM state. 


FIG. 7. (Color online) Energy difference between collinear 
AFM and checkerboard Neel AFM states varying with the 
positive (p) or negative (n) zero-point atomic displacements 
of vertical phonon modes for monolayer FeSe on SrTi 03 in (a) 
undoped and (b) 0.2-electron doped cases. Grey areas high¬ 
light the regions with energy difference within ±10.0 meV/Fe. 
The phonon modes calculated at the equilibrium structure of 
the checkerboard AFM Neel state. 


monolayer FeSe on SrTiOs induced, respectively, by the 
zero-point atomic displacements of all the independent 
vertical phonon modes. The phonon modes are calcu¬ 
lated at the equilibrium structure in the collinear AFM 
order (Fig. 0 ]) and the checkerboard AFM Neel order 
(Fig. EJ), respectively. Only small even negligible vari¬ 
ations are found (< 0.2 /i#), regardless of whether the 
FeSe layer is in a collinear AFM state or in a checker¬ 
board AFM Neel state, nor does it matter what the dop¬ 
ing level of the electrons is. Here the mode numbered by 
32 is exactly the same vibrational mode relating to the Se 
height from the Fe-Fe plane as found in crystal /3 -FeSe in 
our previous work. 36 In comparison, the variation due to 
the corresponding phonon mode can be as large as 1.2 fis 
in the case of bulk FeSe under pressured 

Figures [6] and 0 show how the calculated energy differ¬ 
ence between collinear AFM and checkerboard AFM Neel 
states varies with the zero-point atomic displacement of 
a vertical phonon mode. The phonon modes are calcu¬ 
lated at the equilibrium structure in the collinear AFM 
order (Fig. [6]) and the checkerboard AFM Neel order 
(Fig. 0, respectively. For the undoped FeSe monolayer 
[Figs. H a ) and Ha)], we see that the energy difference 
remains nearly unchanged (^ -23 meV/Fe) for all the in¬ 
dependent vertical phonon modes. In contrast, when 0.2 
electrons are doped, dramatic changes occur [Figs. Mjo) 
and[7|[b)]. There are five modes that reduce the energy 
difference to a value within ±10.0 meV/Fe (represented 
by the dots in the grey area). As shown in Fig. [ljc) and 


Hd), the transition between the collinear AFM state and 
the checkerboard AFM Neel state can be viewed as local 
magnetic moment transfer between Fe atoms or spin flip 
on Fe. Actually, it is well known that it gives rise to mag¬ 
netic frustration when the collinear AFM and checker¬ 
board AFM Neel states are degenerate energetically. 

The atomic displacement patterns of the above six vi¬ 
brational modes, among which one induces the largest 
local moment variation [mode 32 in Figs. HJb) and[5][b)] 
and five substantially reduce the energy difference be¬ 
tween collinear AFM and checkerboard AFM Neel states 
[dotted in gray region of Figs. [6][b) and [Tfb)], are 
schematically shown in Figure [8] The first one is the 
coherent vertical vibrations of Se bridge atoms. The sec¬ 
ond mode involves vertical vibrations of Fe atoms and 
lateral vibrations of Se atoms. The following four vi¬ 
brational modes are contributed by the Ti and O atoms 
at the interface of FeSe epitaxial film and SrTiOs sub¬ 
strate. The calculated frequency of the vibrational mode 
in Fig. [8][d) is 99 meV, which is very close to an en¬ 
ergy difference of 100 meV between two replica electronic 
bands observed in the ARPES experiment^ The vibra¬ 
tional mode with a frequency of 53 meV shown in Fig. 
He) is the ferroelectric phonon. 23,30 It is the vibrational 
mode that could enhance the energy scale of Cooper pair¬ 
ing and even change the pairing symmetry, proposed by 
theorists^ We note that the frequency 53 meV of the 
ferroelectric phonon mode is almost one-half of the fre¬ 
quency 99 meV of the vibrational mode in Fig. H^) and 
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FIG. 8. (Color online) Atomic displacement patterns for six 
typical phonon modes [one mode numbered 32 in Figs. |U[b) 
and[5][b) while five modes dotted in grey region of Figs. [6jb) 
and [?l[b)] of monolayer FeSe on SrTi 03 in 0.2-electron doped 
case. The arrows represent the direction and amplitude of 
atomic displacements. 


meanwhile is close to the energy distance 50 meV of a 
weaker replica band from the original band also observed 
in the ARPES experiment^ In the electron-doped case, 
these vibrational modes effectively alter the interactions 
between Fe atoms directly [Fig. [5fa)-(c)] or indirectly 
[Fig. [8][d)-(f)] and affect the magnetic properties of FeSe 
epitaxial film. 


B. Bilayer FeSe film on SrTiOs 

In the STM experiment, it was found that the mono- 
layer FeSe epitaxial film shows a signal of superconduct¬ 
ing transition around 77 K while the bilayer FeSe epitax¬ 
ial film does not show a superconducting gap. 7 It is thus 
strongly desirable to find out the difference between the 
monolayer and bilayer FeSe epitaxial films. Our previ¬ 
ous calculations indicate that the electronic band struc¬ 
tures of monolayer and bilayer FeSe on SrTiOs in the 
undoped case are very similar^ In the present study, 
we devote ourselves to the effects of both electron dop¬ 
ing and phonon vibrations on the magnetic properties of 
FeSe epitaxial films. 

Figure [9] shows the energy differences between differ¬ 
ent magnetic states for bilayer FeSe on SrTiOs with and 
without electron doping. For the bilayer FeSe, we have 
considered four different magnetic orders, namely i) both 
layers in nonmagnetic state, ii) both layers in collinear 
AFM order, iii) both layers in checkerboard AFM Neel 
order, and iv) the bottom layer in checkerboard AFM 



FIG. 9. (Color online) Energy differences between differ¬ 
ent magnetic states for bilayer FeSe on Ti02-terminated 
SrTiO3(001) surface as functions of electron doping. 


Neel order while the top layer in collinear AFM order. 
The energy of both layers in checkerboard AFM Neel 
order is always energetically lower than that of the non¬ 
magnetic state. The energy difference between both lay¬ 
ers in collinear AFM and both layers in checkerboard 
AFM Neel orders gradually decreases when electrons are 
doped. However, unlike in the monolayer case (Fig. [2]), 
the 0.2-electron doping does not change this energy dif¬ 
ference substantially in the bilayer case. Actually the 
energy difference does not decrease to -17 meV/Fe until 
0.4 electrons are doped. Overall, comparing Fig. [9] with 
Fig. [2j we see that twice as many electrons are required 
to be doped in bilayer FeSe to induce the same effect as in 
monolayer FeSe. This agrees with the recent ARPES ex¬ 
periment, which showed the bilayer FeSe film on SrTiOs 
is much more difficult to be doped with annealing than 
the monolayer FeSe *22 Furthermore, we would like to ad¬ 
dress another important difference between the bilayer 
and monolayer FeSe epitaxial films, namely that the bi¬ 
layer FeSe can adopt different magnetic orders for the 
bottom and top layers. When the electron doping in¬ 
creases (> 0.35e), the energetically favorable state trans¬ 
forms from both layers in collinear AFM order into the 
state of the bottom layer in checkerboard AFM Neel or¬ 
der and the top layer in collinear AFM order. This means 
that enough electron doping alone can induce magnetic 
frustration in the interface FeSe layer of multilayer FeSe 
epitaxial films. 

In Figure [lOj we plot the calculated charge difference 
density for bilayer FeSe on SrTiOs. The electron accumu¬ 
lation at the interface is similar to that in the monolayer 
FeSe epitaxial case (Fig. [3j), meanwhile there is no elec¬ 
tron accumulation in the region between the two FeSe 
layers. Around Fe atoms, there are also charge redistri¬ 
butions. As reported in Tables I, when more than 0.2 
electrons are doped in the bilayer FeSe case, the distance 
between the bottom FeSe layer and the termination Ti02 
layer of the substrate changes little. This facilitates the 
electron doping for the bottom layer in the bilayer FeSe 
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FIG. 10. (Color online) (a) Three dimensional and (b) one di¬ 
mensional charge difference density for bilayer FeSe adsorbed 
on SrTi 03 with 0.2-electron doping. Dark red and blue iso¬ 
surfaces in panel (a) are respectively electron accumulation 
and depletion areas in isovalue of 0.01 e/A 3 . 


case. On the other hand, the distance between the bot¬ 
tom and top FeSe layers tends to increase when the elec¬ 
tron doping exceeds 0.2 e, which makes the top layer hard 
to be further doped. For the monolayer and bilayer FeSe 
films, the similar electron accumulation at the interface 
suggests that the similar physical phenomena would oc¬ 
cur at the inter facial FeSe layer. 


IV. DISCUSSION 

In iron-based superconductors, it is commonly 
thought that AFM fluctuations are responsible for the 
superconductivit y 4Q i 41 For bulk [3 -FeSe under pressure, 
the zero-point atomic displacement of the vertical vibra¬ 
tional mode of Se atoms induces a large variation on lo¬ 
cal magnetic moment via strong spin-phonon coupling, 
as demonstrated in our previous work. 36 - For the FeSe 
epitaxial films in the present study, it is found that, un¬ 
like bulk FeSe^ the variation of local magnetic moment 
via zero-point atomic displacements does not show any 
meaningful enhancement in both undoped and doped 
cases (Figs. 0] and [5j). In contrast, the energy differ¬ 
ence between collinear AFM and checkerboard AFM Neel 
states shows a meaningful reduction with electron doping 
for both the monolayer (Fig. [2]) and bilayer FeSe (Fig. 
[9j) cases. 

For monolayer FeSe on SrTiOs, the dramatic reduc¬ 
tion in the energy difference between the collinear AFM 
and checkerboard AFM Neel states is caused by a com¬ 
bined effect of electron doping and phonon vibrations. 
As shown in Fig. [2j without considering the phonon ef¬ 
fect, the energy difference decreases with electron dop¬ 
ing but saturates at about -17 meV/Fe after more than 
0.2 electrons are doped. On the other hand, without 


including electron doping, the absolute value of the min¬ 
imum energy difference is still larger than 23 meV/Fe 
when only the phonon effect is taken into account [Figs. 
Ha) and Ha)]. The combined effect of electron doping 
and phonon vibrations brings about a certain probability 
to substantially reduce the energy difference [Figs. 
andHb)]- I n the J 1 -J 2 Heisenberg model^ the mag¬ 
netic coupling J\ between the nearest-neighboring (NN) 
Fe spins and the J 2 between the next-nearest-neighboring 
(NNN) Fe spins can be derived from the energy differ¬ 
ences between different magnetic states. In the undoped 
case, it gives 29.0 meV/5' 2 for J\ and 20.8 m eV/S 2 for 
J 2 with S being the local magnetic moment on Fe. With 
electron doping and phonon effects, it is very likely that 
the reduced energy difference between collinear AFM 
and checkerboard AFM Neel states makes J 2 compara¬ 
ble with Ji/2. Then the magnetic frustration or mag¬ 
netic instability takes place, which is helpful to super¬ 
conductivity. This provides a consistent physical picture 
for previous experimental^ - — and theoretical studies on 
the FeSe monolayer^ - — 

The bilayer FeSe epitaxial film behaves distinctly 
from the monolayer one. At regular doping levels, the 
Fermi surfaces of monolayer and bilayer FeSe detected in 
ARPES experiments are quite different^ - - 3 The former 
shows only electron-type Fermi surfaces around the Bril- 
louin zone corners, while the latter shows both electron- 
type Fermi surfaces around the corners and hole-type 
Fermi surfaces around the center of the Brillouin zone. 
Nevertheless, when it is effectively doped with prolonged 
annealing time, the bilayer FeSe epitaxial film demon¬ 
strates a similar Fermi surface to the monolayer one in 
ARPES measurement. 20 This means that the top layer 
of bilayer FeSe is much more difficult to be doped than 
its bottom layer experimentally. On the other hand, we 
remind the reader that the Fermi surface of monolayer 
FeSe observed in ARPES can be well reproduced in cal¬ 
culations when the FeSe layer is set in checkerboard AFM 
Neel order j 24 i 25 Meanwhile, our calculations show that to 
make the checkerboard AFM Neel order favorable ener¬ 
getically by doping, the top layer of bilayer FeSe is much 
harder than its bottom layer (Fig. [9|). This is consistent 
with the ARPES experiment. 20 In addition, in transport 
experiments, clear superconducting transitions have been 
observed both on a five unit-cell FeSe film covered with 
an amorphous Si protection layer— and on a monolayer 
FeSe film protected with FeTe film.— It has been argued 
that the superconductivity occurs at the first unit-cell 
FeSe in transport measurements^ 15 In the calculations, 
we have found that the magnetic frustration in the inter¬ 
facial FeSe layer of bilayer FeSe film may be induced by 
electron doping, which supports this viewpoint. 


V. SUMMARY 

We have studied the magnetic properties of mono- 
layer and bilayer FeSe epitaxial films on SrTiOs with 
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electron doping by using the first-principles calculations 
with vdW correction. For monolayer FeSe epitaxial film, 
the combined effect of electron doping and phonon vibra¬ 
tions gives rise to a certain probability to induce magnetic 
frustration between the collinear AFM and checkerboard 
AFM Neel states, which is helpful to superconductivity. 
Most of the effective phonon modes, which can substan¬ 
tially reduce the energy difference between AFM states 
with electron doping, are contributed by interfacial Ti 
and O atoms. In contrast, for bilayer FeSe epitaxial 
film, the interfacial FeSe layer has magnetic frustration 
induced by electron doping much more easily than the 
top layer, which suggests that for multilayer FeSe epi¬ 
taxial films, the superconductivity readily takes place at 
the bottom layer. These calculated results agree with the 
existing experimental studies. 
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Appendix 

The relaxed structural parameters for monolayer and 
bilayer FeSe epitaxial films in the nonmagnetic state are 
listed in Table II. 


TABLE II. The relaxed layer distances d (in A) of monolayer 
and bilayer FeSe epitaxial films in the nonmagnetic state un¬ 
der different levels of electron doping. 


doping(e) 

0.0 

0.1 

0.2 

0.3 

0.4 

Monolayer 

dTiO-FeSe 

4.138 

4.134 

4.137 

4.137 

4.135 

Bilayer 

^TiO-FeSe 

4.154 

4.138 

4.143 

4.140 

4.134 

^FeSe-FeSe 

5.098 

5.127 

5.127 

5.132 

5.146 
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